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ZES Z YT 


i 


K. URBANIK, A. ZIIJBA 


Prognoza akty wnosci slonca 

I RqKopis dostarczono 18. 11. 1955. 

Zmiany elektromagnetycznej akty wnosci slonca, powoduj^ce w szcze- 
golnosci zaburzenia w odbiorze radiowym, charakteryzuje si^ na ogol licz- 
bami Wolf a, ktore opisuj^ ro2anieszczenie i ilosc plam slonecznych. Poniew^ 
scisla zaleiznosc liczb Wolf a od czasu nie jest znana, przyjmuje si^, ze 
tworz^ one proces stochastyczny. W artykule bmowione s^ krotko probabili- 
styczne metody prognozy liczb Wolf a: liniowa metoda Kolmogorowa-Wienera 
oraz metoda uwzgl^dniaj^ca koszt eksperymentu oparta na teorii funkcji 
decyzyjnych. 

. Akty wnosci^ elektromagnetyc 23 n^ slonca interesuj^ si^ zarowno astro- 
nomowie i fizycy jak i radiotechnicy. Jej znaczenie dla radiotechniki po- 
lega na tym, ze w duzym stopniu od niej zalezy dzialariie shizb telekomu™ 
Vakacyjnych. Zwi^kszenie aktywnosci wplywa ujemnie na jakosc odbioru. 

Aktywnosc slonca jest scisie zwi^zana z wystgpowaniem plam slonecz- 
nych. Dlatego opisujemy za pomoc^ tak zwanych liczb Wolfa, charak- 
teryzuj^cych ilosc i rozmieszczenie plam na slohcu 

r^k (lOp-l /) 

- 1 - 

r — jest tu liczbq Wolfa, 

g — oznacza ilosc grup plam slonecznych (izolowane plamy uwazamy 
za oddzielne grupy), 
f — calkowit^ ilosc plam; 

k — jest wspolczynnikiem zaleznym od przyrz^du pomiarowego,. 

w przyblizeniu rownym jednosci..^.. 

Zaleznosc liczb Wolfa od czasd nievjest dokladnie znana. 

Obserwacje wskazuji na okresowe nasilanie si^ liczby plam slonecz- 
nych, przy czym wyrazny jest tu cykl jedenastoletni. 

W pierwszym przyblizeniu przyjmuje si^, ze liczby Wolfa zalezy od 
czasu ^losowo. Pozwala to na wykorzystanie aparatu teorii procesow sto- 
chastycznych do badania aktywnosci slonca, a w szczegolnosci do ustalania 
prpgnozy aktywnosci slonca. ’ 
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Niech r(t) b^d^ca funkcj^ losow^ oznacza liczb^ Wolfa odpowiadaj^c^ 
chwili t. Obok funkcji r(t) wygodnie jest poslugiwac si^ odchyleniem 0 (t) 
wielkosci r(t) od sredniej wartosci r (t) ‘ 

g(t) = r(t)—r(t) . 

Dla odchylenia e (t), ktore rowniez jest funkcji losowg, mamy 

g(t)=o. 

Przyjmujemy takze, zgodnie z doswiadczeniem, ze korelacja odchylen 
0 (t) i e (^ + zalezy wylacznie od diugosci przedzialu czasowego r, a wi^c 
funkcja korelacyjna ma postac 

I 

Br(T) = Q(t-hr)^(T), 

Odchylenia g (t) tworz^ zatem stacjonarny proces stochastyczny. Przy- 
pusciTiy, ze .znamy liczby Wolfa w chwilach Wowczas szacu- 

j^c r(t) sredni^ arytmetyczn^ 


- YriU) 

jLJ , 


mozna obliczyc odchylenia q (tj) , . . . , ^ (tn) . Zagadnienie prognozy liczby 
Wolfa w chwili t, przy znanych liczbach Wolfa wchwilach tn 

sprbwadza si^ wi^c do prognozy odchylenia Q{t), gdy znane odchylenia 
... Ograniczymy siq do prognozy polegaj^cej na oszaco- 

waniu nieznanej wartosci q (t) kombinacjq liniow^ wartosci ^ (tj) , . . . , ^ (tn) . 
Poza funkcja korelacyjna w prognozie nie wykorzystujemy innych cha- 
rakterystyk procesu stochastycznego git) . 

Teoria prognozy liniowej opracowana zostala przez A. N. Kolmogo- 
rowa, N. Wienera i innych. 

Prognoza liniowa sprowadza si^ do znalezienia ukladu wspolczynnikow 
liczbowych , . . . , an ; odchylenie q (t) szacuje si^ wowczas kombinacj^ 
liniowq (ti) + (tn)* Prognozy opisan^ wspolczynnikami oi , . . , an 

hazywamy optymaln^, jesli sredni bl^d kwadratowy prognozy jest naj- 
mniejszy, to znaczy 


(7^= g(t) — J^ ttjQ (tj) I = min i o (t) - "iO.it) ■ 

i = l ' I . y, , . V . , •!,. i i = l 

- r;'". 

Dowodzi siQ, ze wspolczynniki ai , . . . , a,i opisuj^ce najlepsz^' prognozy 
spelniaj^ uklad rownah (por. [2]) ' ‘ v , 

■ ' ■ - ■ 

^ QiB (ti tfc) = jB (t ■ " tk) (fc = 1 , 2 , . . . , Tl) ■ . 

‘ Kreska nad symbolem zmiennej losowej oznacza wartosc sredni^ (oczSkiwan^), 
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, a sredni bl^d kwadratowy (j^ - optymalnej prognozy jest rowny 

nn 

^ OfcOiBfti— t(c). 

. k=l i=l 

Odchylenia q (tj) (tn) mozemy, jak powiedzielismy, obliczyd. Po- 

zwala to oczywiscie szacowac funkcj^ korelacyjn^ B (t) wedlug znanego 
wzoru (por. [2] lub [6]) 

* » B{r)^ eiti + T) Q{ti) . ■ 

tn ti~ T 

' .r . t, + T<t„ 

Znaj^c zas funkcjQ B(t) z rownan (1) mozemy wyznaczyc wspolczynniki 
optymalnej prognozy. 

Zanim podamy pewne konkretne przyklady prognozy, zwrocimy uwag^, 
ze omawiana .teoria prognozy liniowej podaje rachunkowy sposob otrzy- 




Wbi ' ' f952 

Rys. 1. Wykres aktywnosci slonca wedlug prognozy i na podstawie 
obserwacji w latach 1951/1952. 


mania optymalnej prognozy (w podanym wyzej sensie) na podstawie 
d a n y c h obserwacyjnych. Material obserwacyjny aktywnosci slonca 
jest niezmiernie bogaty i nie moze bye caly wykorzystywany w rachun- 
kach liczbowych. Powstaje zatem zagadnienie wyboru optymalnej i 1 o- 
s c i obserwacji i ich optymalnego rozmieszczenia. Innymi slowy 
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powstaje pytanie — ile i ktore obserwacje nalezy wykorzystac. Za- 
gadnienie to jest nierozwiqzane. Proba, podj^ta przez autorow tej pracy, 
wskazuje, ze zagadnienie prowadzi do rownan, ktorych nawet bardzo nie- 
dokladne rozwigzanie jest praktycznie n'a razie niemozliwe. (Do zagad- 
nienia tego powrocimy jeszcze przy omawianiu innych teorii prognozy). 
W teorii prognozy liniowej przyjmuje si^ heurystycznie takze (por. [3], 
[4]) wzcry typu . .. 

e(t)=aie(t-l)+a2e(t-2) + a,oe(t-lO)+aii0(t-ll) + 

'V 

+ 020? (t— 20) + a2i5 (tN;-21)+. 

Na przyklad, dla dwoch obserwacji otrzymuje si^ na r(t) wzor 
r (t)= 14,5 + 1,3 r (t- 1)- 0,6 r(t- 2). 

Prognozy na lata 1925—1929 opartq na obserwacjach z dwoch poprze- 
dnich lat przedstawia ponizsza tabhca (por. [3]). 


Rok 

Prognoza 

Obserwacja 

Bt^d 

1925 

33,2 

44,3 

11, "i 

1926 

63,3 

6^9 

0,6 

1927 

71,5 

69,0 

-2,5 ^ ’ 

1928 

65 ;4 

77,8 

12,4 

1929 

^ ^ 74,0 

65,0 

—9,0 


Uwzgl^dnijmy teraz tzw. koszt wykorzystania danych obserwacyj- 
nych C(N). Koszt ten nalezy rozumiec bardzo ogolnie; wchodzi do niego 
zarowno czas potrzebny do obliczeh jak i zuzycie maszyn do liczenia itp. 
Oczywiscie zalezy na zmniejszeniu tego kosztu. Niech Qn oznacza prognozy 
odchylenia ? (t) opart^ na N obserwacjach (obserwacje te traktujemy jako 
zrnienne losowe, a wi^c gw jest wielkosci^ losow^). W(? ,?jv) oznacza do- 
kladnosc prognozy odychlenia liczby Wolfa od sredniej wartosci. Moze tp- 
byc na przyklad 

w(?,?jv)=!?— ? n|- 

Zalezy rowniez na zmniejszeniu funkcji W (?, ?Ar),to jest na polepszeniu 
dokladnosci prognozy. 

Gdybysmy umieli oszacowac strata gospodarcz^ wynikaj^c^ z danej 
bl^dnej prognozy, wowczas mozna by otrzyma6 nowe wyniki w zakresie 
poruszonego juz zagadnienia ilosci i rozrhieszczenia w czasie danych ob- 
serwacyjnych. Dokladne oszacowanie tej straty jest przypuszczalnie tru- 
dne. Jest to oczywiscie trudnosc pozamatematyczna. -Pohiewaz jednak 
w praktyce wystarczy oszacowanie przyblizone, podamy przyklad wyzna- 
czania ilosci obserwacji w przypadku, gdy funk'cja C(W) i W(?,?iv)s^ zna- 


.vj 
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nymi funkcjami, ktorych wartosci wyrazaj^ si^ w zlotych.' Wskazemy 
zreszt^ pozniej pewne warianty tej teorii oparte na innych (w pewnych^ 
przypadkach dogodniejszych) zalozeniach o funkcji W(g , Pn) • 

Przejdziemy teraz do teorii majqcych na celu prognozQ dystrybuanty 
procesu. Rozpoczniemy od w^skiego zagadnienia szacowania ilosci obser- 
wacji. Pozniej zajmiemy si«j zagadnieniami 'o wiele ogolniejszymi. 

Sume 

Rq(.N ,qn)=C {N\+W {q , qn) 

nazywamy calkowitym kosztem prognozy, ^a wyrazenie 

R (iV , eN) = maxiC (iV)+W (e , ew)] 

' / ® ■ . ■ 

maksymalnym calkowitym kosztem prognozy. Hose obserwacji potrzeb- 
nych do prognozy nazywamy optymaln^, jezeli 

R(No,QN)^mmRiN ,on). 

■N : ■ 

Na prostym, szczegolnym przykladzie pokazemy metod^ znajdywania 
optymalnej ilosci obserwacji. 

Przyjmiemy najpierw pewne zalozenia upraszczajqce. Przede wszyst- 
kim zakladamy, ze obserwaeje wykonujeniy w sposob ci^gly, a nie dys- 
kretny. Dysponujemy wiec wykresem rejestrowanym przez aparat od 
pewnej chwili pocz^tkowej 0. Dalej zakladamy, ze odchylenie q (t) liczb 
Wolfa od sredniej tworz^ proe'es normalny, czyli ze rozklady prawdopo- 
dobienstw sq dane wzorem 


P[e(t)<a:] = 


i-fi 

2'KtaJ i 




i ze prognoza dotyezy tylko niewiadomego parametru /« w tym rozkladzie. 
Przyjmijmy, ze dokladnosc W (/<, A) prognozy, jest funkeja postaci .. 

W(,«, A) = to(.| /t-A|), 0 

gdzie ! . . ' -'tf-'r-H 

«— jest prawdziw^, ‘ ■ v. ■ 

A— przepowiedzian^ wartosci^ parametru. ' . * 

Prognoza wielkosci [jl polega na oszacowaniu.za pomoc^ sredniej arytme- 
tycznej wynikow. obserwacyjnych. c‘ . \ 

Dojwodzi siQ [1], ze optymalna ilosc T pomiarow przy pe^nych zalo- 
zeniach regularnosci funkcji w{x) i C(x) (ktore w praktyce mozna przyj^c) 
otrzymuje siQ z rozwi^zania rownania “ ” " ’ !. 


dC (T) 


w(x) e 
\/2 tiT • a 


w (x) x^ 
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wzgl^dem T. Znaj^c funkcje C(T) i u>(x) mozna to rownanie rozwi^zac » 

metodami przyblizonymi. 

Podamy teraz pewien inny typ prognozy, maj^cy o wiele szerszy za- 
si^g niz podana na przykladzie teoria wyznaczania optymalnej obser- ' 

wacji. Teoria ta, oprocz optymalnej ilosci i optymalnego rozmieszczenia * 

pomiarow, pozwala na skonstruowanie optymalnej prognozy [5]. 

Podstawowym poj^ciem tej teorii jest poj^cie funkcji decyzyjnej. 

Przez funkcjQ decyzyjnq rozumiemy ci^ przepisow, z ktorych pierwszy 
orzeka, jakie obserwacje powinnismy wykonac i uwzgl^dnic na pocz^tku, '* 

drugi — na podstawie tych wynikow — b^dz orzeka, ze wykonujemy 
i uwzglQdniamy dalsze obserwacje (i jesli tak, to jakie), b^dz wydaje 
tzw. decyzje ostateczn^. Trzeci przepis wchodzi w rachub^ tylko wtedy, 
jesli drugi nie wydal decyzji ostatecznej i znowu, na podstawie wynikow 
obserwacji dokonanych na obu poprzednich etapach b^dz zaleca dalsze 
pomiary, b^dz wydaje decyzje ostateczn^ itd. Przez decyzje ostateczn^ 
rozumiemy orzeczenie, ze badany proces stochastyczny ma tak^ a tak^ 
dystrybuant^. W samym poj^ciu funkcji decyzyjnej abstrahujemy od 
tego czy jej zalecenia s^ dobre, czy decyzja ostateczna jest prawdziwa 
itp. Przez dystrybuante procesu rozumiemy funkcje 

I » • • • > tn (®i > • • • » ^n) P {Q (ti) ^ 1 • ' • } Q i^n) ^ Xn} 

okreslan^ dla wszystkich uktadow ,ti , . . . , t„ , Xj , . . , , x« , gdzie symbol 
P{ } oznacza, jak zwykle, prawdopodobienstwo zdarzenia opisanego w na- 
wiasie. 

Uzywane w teorii funkcji.. decyzyjnych badanie etapowe nazywamy 
postepowaniem sekwencyjnym. , 

Oprocz wyzej opisanych funkcji decyzyjnych uzywamy funkcji decy- 
zyjnych zrandomizowanych (losbwych) [5], Ma to zastosowanie do przy- , 
padkow, gdy prognozy konstruujemy nie raz, lecz wiele razy. Mozemy ^ 

wtedy w poszczegolnych badaniach stosowac rozne funkcje decyzyjne. 

Jesli, ogolnie, z danego ukladu funkcji decyzyjnych za pomoc^ ustalo- 
nego mechanizmu (kostki) bedziemy losowali poszczegolne funkcje i sto- ■ 

sowali je do kolejnych prognoz, b^dziemy mowili, ze posiugujemy sie \ 

funkcjq decyzyjne zrandomizowane. 

Zalety randomizacji mozna stwierdzic na prostym, nawet nie sekwen- 
cyjnyni przykladzie. Ktos kladzie monete, my dajemy prognozy poja- 
wienia sie orla. Z’adanie powtarza sie wiele razy. Jesli b^dziemy stale po- 
dejmowali jedne i te same decyzje z dwu mozliwych: pojawi 
sie orzel, pojawi sie reszka, to, jezeli moneta byla rzucana przypadkowo, 
mozemy gwarantowac, ze srednio co najmniej (a takze co naj wyzej) po- 
lowa prognoz bedzie dobra. Jesli posluzymy sie decyzje zrandomizowane 
i sami swoje odpowiedzi: orzel, reszka, bedziemy losowali monete, to 
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takze b^dziemy mogli zagwarantowac, ze srednio poiowa naszych prognoz 
b^dzie prawdziwa, ale juz nie bedziemy potrzebowali za- 
kladac ,jak poprzednio, ze moneta, ktorej rzuty mamy 
odgadywac, jest rzucona przypadkowo. Nie bedziemy potrze- 
bowali w ogole zadnych hipotez na ten temat. Randomizacja ma jeszcze 
inne zalety, ktorych nie bedziemy tu szczegolowo rozwazali. 

Przypuscmy, ze mamy dan^ rodzin^ kompletnych dystrybuant pro- 
cesow stochastycznych. Dystrybuanty oznaczymy symbolami Fi ^ Fn- 
Tq dan^ rodzin^ ' bedziemy nazywali baz^. Utworzmy zrandomizowanq 
kombinacjQ liniow^ 

cctiFti (a^ ^ 0 , ^Cfj= 1) (2) 

bazy. Dzi^ki normalizacji, tak^ ' kombinacj^ znow mozemy uwazac za 
dystrybuanty. 

Przypuscmy, ze tak wybralismy bazy dystrybuant (w praktyce kilka- 
dziesi^t), iz mamy prawo .przypuszczac, ze losowy proces aktywnosci 
slonca daje siy w przyblizeniu opisac jak^s z dystrybuant postaci (2). 
Przypuscmy dalej, ze wybralismy uklad funkcji* decyzyjnych di , . . , , dm 
(w praktyce — kilkadziesi^t sekwencyjnych metod testowania hipotez, 
ze proces q (t) ma dystrybuanty F{, F2 > • • * > * Na ich podstawie skon- 

struujemy optymalna funkcjy decyzyjn^ do badania procesow ksztaltu (2). 

Niech C{N) oznacza, jak poprzednio, koszt wykorzystania N obserwacji, 
W(FiyFj) — straty gospodarcz^ wynikaj^c^ z powziycia decyzji ostatecz- 
nej orzekaj^cej, ze proces Q(t) opisany jest dystrybuanty Fj, podczas gdy 
naprawdy opisany jest dystrybuanty Ft [W (Ft , F^) > 0 , W {Ft , Ft) = 0] . Przy- 
puscniy na chwily, ze proces aktywnosci slohca gt opisany jest dystry- 
buanty F^ i ze poslugujemy siy funkcjy decyzyjny Znajyc dystrybuanty ^ 
Fi i funkcjy decyzyjny d^ la two mozemy obiiczyc sredni koszt calkowity ^ 
badania za pomocy funkcji d^, procesu opisanego dystrybuanty Ft , koszt 
ten wyrazi siy wtedy funkcjy 

R {Fi , djc) . 

Zadaniem naszym jest ten koszt uczynic jak najmniejszym. 

Jesli proces aktywnosci slohca opisuje siy dystrybuanty 

• «iFi+ . . . . , -\~CtnFn , 

a my poslugujemy siy funkcjy zrandomizowany 

/^l^l , . . . , ^mdm (3l 

{^i oznacza prawdopodobiehstwo wybrania funkcji decyzyjnej di) , to sre- 
dnio koszt calkowity wyrazi siy formy biliniowy 

• “ Koszt wykorzystywania pomiarow i straty wynikajqc^ ze zlej decyzji. 
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i=n , j = m • 

. ^ d,), 

^=1 » j = i 

ktor^ dla skrotu oznaczymy ' 

Rijf'Hpj . 

Maksimum formy Rijai^j ze wzgl^du na at przy ustalonych Pj nazywamy 
kosztem gwarantowanym przez funkcjQ. decyzyjn^ (3)' (zrandomizowan^). 

FunkcjQ decyzyjn^ zrandomizowanq b^di+ bndn nazywamy 

najlepsz^ jezeli gwarantuje minimalny koszt 

Rljdi^ j . 


max 




Rijaibj = min 


max 


Aby zilustrowac rachunkowe trudnosci przy efektywnym znajdowaniu 
wspolczynnikow b,*, rozwi^zemy prosty przyklad przyjmuj^c i='2, j^2, 
Przy wiekszych wartosciach i, j trudnosci wzrastaja podobnie jak przy 
przechbdzeniu od dwoch rownan liniowych z dwiema niewiadomymi do 
wi^kszej liczby rownan i niewiadomych. 

Frzyjmijmy RijajPi = ajP^ + a2p.^, ' ' ^ 

Mamy ' wyznaczy c min max + 

PljPi « 1,«2 

Ekstremum mieszane minimum maximorum wyznaczamy podobnie jak 
zwyczajne, przez przyrownanie pochodnych' do zera. Aby ulatwic sobie 
analityczne operowanie zmiennymi ktore s^, jak wiemy do- 

datnio unormowane najlepiej zamiast tych zmiennych uzyc zmiennych 
;sin^a, cos‘^a , sin^/? , cos^p , 

Rozniczkujac form^^-^ 

sin^a • sin^jS H-cos'-^a • cos"/5 , 
otrzymujemy na ekstrema uktad rownan 


(4) 


2 {sin^P — cos'^P) 
2 (sin^a — cos^a) 


sin a cos a==0, 
sin P cos p = 0 . 


TJklad ten ma pierwiastki: 


1) sin a— sin ^ = 0 , 

2) cos a = cos = 0 , 

3) sin a = cos ^ = 0 , 

4) cos P^sin a = 0, 

d) sin^a = cos^a=sin^P = cos^P = — . 

2 

Latwo sprawdzic, ze pierwiastkom 1) i 2) odpowiada maksimum, pier- 


To znaczy 


ai + a2==l, cii ^ 0 , 02^0, 

+ = ft>0, ft>0. 
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^ j ^ 

K wiastkom 3) i 4) — minimum formy 4. Pierwiastkom 5) odpowiada ekstre- 

, . , 1 , 1 

mum mieszane. Znalezlismy wi^c ekstremalne wartosci Oi — , Rz— . 

ij ■ ■ . , 

' Wynik ten mozna bylo przewidziec, gdybysmy przedtem wyjasnili, ze 

forma + odpowiada zagadnieniu prognozy orla i reszki w przy- 

kladzie, o ktorym wspominalismy. Wspolczynniki bi , bj maj^, znaczenie 

prawdopodobienstw, z jakimi powinnismy przepowiadac orla i reszk^. 

Oczywiscie teoria decyzyjna zostala przedstawiona w sposob bardzo 

pobiezny i szkicowy. ' 

'Wspomnijmy jeszcze, ze nast^puj^ce drogi orniniQcia trudnosci sza- 
cowania straty gospodarczej ze zlej prognozy. Koszt wykorzystania ba- 
dan i jakosc prognozy mozna mierzyc w rdznych jednostkach i powyzszq 
teori^ iiprawiac jak gdyby w zakresie liczb zespolonych. IVtedy zamiast 
optymalnej funkcji decyzyjnej otrzymujemy rodzin^ optymalnych nie- 
porownywalnych mi^dzy sob^ funkcji, ostatecznego zas wyboru doko- 
nujemy w zaleznosci od lokalnych potrzeb i mozliwosci. Teoria dostarcza 
jak gdyby ofert, o ktorych wie, ze sq lepsze od innych mozliwosci, ale 
mi^dzy ktorymi juz nie rozstrzyga. 

, '' Innq drog^ jest budowa analogicznej, cHociaz istotnie roznej teorii 

opartej, z grubsza mowiqc, na tnozliwosciach finansowych instytucji pro- 
wadz^cej badania. Chodzi o to, ze w praktyce i tak na takie a takie bada- 
f dania przeznaczamy tyle a tyle. Teoria wskazuje natomiast optymalny 
schemat badaA opartych na danym budzecie i, co wiQcej, pozwala sekwen- 
cyjnie badac czy dalsze zwiQkszenie budzetu jest celowe, czy tez budzet 
ten nalezy zmniejszac. - 

Trzeba zaznaczyc, ze przy stosowaniu wszystkich wyzej przedstawio- 
nych teorii prognoz w praktyce napotykamy dose duze trudnosci. Wynika 
to z faktu, ze chociaz teorie te s^‘ na ogol szeroko i precyzyjnie rozbudo- 
wane, to jednak prowadz^ czQsto do niezmiernie skomplikowanych- ra- 
chunkow liumerycznych lub wymagaj^ trudnych i kosztownych badan. 

Instytut Matematyezny 
Pnlahiej Akademii Nauk 
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K. yPBflHHK, a. 3EM5a 

, ' nporHOs cojihehhoh aktmbhoctm 

\ 

P e 3 K) M e / 

MsMeneHHH ajieKTpoMarHMTHOw cojinenHOM aKTMBHocTw, BbiabieaioiUMe b ocoben- 
HOCTM noMexM npMeMa paAMonepe^anw, xapaKTepwayiOTCH Boobiqe HMCJiaMM . Bojibcba, 
oriMCbiBaioinMMM pacnpe/^ejieHwe pi KOJiMnecTBo cojineHHbix nATen. IIoaTOMy npornos 
COJIHeHMOM aKTPIBHOCTM CBOAMTCH K nporH03y HMCCJI BOJIbC^a. TOHHaH 3aBPlCHMOCTb 
HPiceji Bojibcba ox BpeivieHPi — naM HepiSBecxHa. B nepsoM nppi6jiPiJKeHHPi nppiHM- 
Maexcfl, HTO MSMeneHHH nnceji Bojibc^a mmciot cjiynanbiM xapaKxep. 3 to nosBOJinex 
picnojibaoBaxb B npornose annapax cxoxacxMHecKMx ^po^eceoB. B cxaxbe BKpaT^e 
o6cyxc;^eHbi npoSabMJiPicxpiHecKPie Mexo^bi npornosa HHCJia Bojibcha. Jlynuie Bcero 
Mccjie/^oBaH xeopexMHecKM Mexo^ jiwHeMHoro npoxHoaa KojiMoropoBa - Bpinepa, koxo- 
pbiM xeM Jiynuie, neu bojibmee hmcjio HaSjiioAaxejibHbix ;^aHHbix oh Mcnojibsyex. 
OAHaKO 6oAbmoe hmcjio nabjiioflaxejibHbix 7^aHHbIx Be^^ex k ajimhhbimm m KponoxjiM- 
BbiM BbiHMCJieHHM. Ho cjiy^yex noaxoMy npcHebperaxb Mexo^aiviM npornosa, ncnojib- 
ByiomMMM Hebojibiiloe hmcjio ^aHHbix m ^aiomMMM cpaBHMxejibHo xopouiMe peayjib- 
xaxbi. B cxaxbe npi^se^eH xbkom Mexo^(, ynMXbiBaiomMM cxommocxb onbixa m ocHOBaH- 
HblH Ha XeopMM ;^eUP13M0HHbIX CjjyHKMMM. 

K. URBANIK, A. ZI^BA 

PREDICTION OF SOLAR ACTIVITY 
® Summary 

Electromagnetic, changes in solar activity, which cause disturbances in radio re- 
ception, are characterized in a general manner by Wolf numbers. These numbers 
describe the location and number of sunspots. Thus, prediction of solar activity is 
reduced to the prediction of Wolf numbers. TJie exact dependency of Wolf numbers 
on time is not known. As a first approximation, it is assumed that the change in the 
Wolf numbers is of a chance nature. This permits stochastic process apparatus to 
be used for predictions. In this paper, probalistic methods for predicting Wolf num- 
bers are briefly discussed. The method most investigated is the Kolmogorov-Wiener 
linear method of prediction. In this method, the greater the amount of data amploy- 
ed, the better the results obtained. A great quantity of observed data, however, lead 
to long and timeconsuming calculations. Therefore, prediction methods employing 
not too great quantitative data, and giving relatively good results, are not without 
practical importance. In this paper, one such method, based on the theory of decision 
functions, is described. The cost of the experiment has also been taken into account. 
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621.396.11 

S. BOROWSKI, S. JASmSKI, S. MANCZARSKI 

Zjawisko Dop^lera w propagacji jonosferycznej 

RqJcopis dostarczono 13. 12. 1955 

Praca podaje zestawienie wynikow pomiarowy'ch odnosnie zjawiska 
Dopplera w oswietleniu klasycznej teorii tego zjawiska. Poniewaz niektore - 
wyniki doswiadczalne zmian cz^stotUwosci przypisywane zjawisku Dopplera 
znacznie odbiegaj^ od obllczeri teoretycznych, przedstawiono przypuszczalny 
mechanizm zmian cz^stotliwosci w warunkach granicznych (bliskich MUF) 
odbic od jonosfery. - 

1. WSTFP 

Celem pracy jest wyjasnienie rozbieznosci poglqdow, jakie spotyka 
si^ w literaturze na temat znaczenia w radiokomunikacji krotkofalowej 
zjawiska Dopplera, ktore zachodzi w zwi^zku z niestabilnoscig jonosfery. 
Cytowane przez roznych aiitorow (Dodatek 1) pomiarowe wyniki wahan 
czQstotliwosci powodowanych przez to zjawisko rozni^ si^ mi^dzy sob^ 
o kilka rz^dow wielkosci. Niektorzy autorzy wi^z^ wystQpowanie silnego 
efektu Dopplera z duz^ zmiennosci^ k^ta przychodzenia fall jonosfe- 
rycznej. 

• Wedlug badan przeprowadzonych na falach ultrakrotkich (ponizej 
150 MHz) w ostatnich latach w rejonach zorzy polarnej stwierdzono, ze 
przy odbiciach od kotary zorzowej wyst^puj^ dose znaezne zmiany czQsto- 
tliwosci, ktore mog^ by6 przypisane zjawisku Dopplera. Podezas jednej 
z serii prob przy 50 MHz otrzymano zmiany cz^stotliwosci fali nosnej 
zawarte w granicach od 0 do 200 Hz [6] . 



2. KLASYCZNA TEORIA ZJAWISKA DOPPLERA 
NA TLE WYNIKOW DOSWIADCZALNYCH 

Zjawisko Dopplera znane w akustyce i w optyce (przy zjawiskach 
astronomicznych) obserwowane jest rowniez przy propagacji fal radiowych 
odbijanych od jonosfery. 

Zjawisko to polega na tym, ze cz^stotliwosc odbierana rozni si^ od na- 
dawanej. Wedlug klasycznej teorii zjawiska Dopplera obserwowana 
zmiana cz^stotliwosci moze wynikac na skutek: 
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1. ruchu zrodta, 

2. ruchu odbiornika, 

3. ruchu srodowiska (w przypadku propagacji jonosferycznej fal radio- 
wych — wskutek ruchu warstwy odbijaj^cej), 

4. kombinacji ruchow wyszczegolnionych w punktach 1. i 3. 

• Zjawisko' Dopplera przy ruchach zwierciadlanej warstwy jonosferycz- 
nej podane zostalo w uj^ciu optyki geometrycznej w Dodatku 2. Z dys- 
kusji wzorow przytoczonych w tym Dodatku wynika, ze zjawisko Dop- 
plera przy propagacji jonosferycznej moze bye wywolane jedynie przez 
skladow^ pionow^ szybkosci poruszania siQ warstwy odbijaj^cej. Zmiana 
cz^stotliwbsci (wzgl^dna) wyraza si^ wowczas wzoreni 

zd/ u . 

— =2 — sma, (1) 

J c 

gdzie 

— zmiana, cz^stotliwosci (na skutek zjawiska Dopplera), 
f — czQstotliwosc emitowana, 

u — skladowa pionowa szybkosci unoszenia si^ warstwy, ■ 
c — szybkosc swiatla, 
a — k^t elewaeji. 

Ze wzoru (1) wynika, ze ^ jest tym wiQksze, im wi^ksza jest szybkosb 
u i k^t a . 

W tablicy 1 podano wyniki obliczeh ~ dla k^tow elewaeji 10° i 30° 

na podstawie wzoru (1) przy obserwowanej, na przyklad w Anglii w okre- 
sach malych plam slonecznych, szybkosci unoszenia si^ warstwy odbija- 
jqcej 12 m/sek. 

T a b 1 i c a 1 





a"" 

u m/sek 

y -lO® 

10 

! ^2, 

1,4 

30 

i 

4 


Z porownania powyzszyeh obliczeh z danymi Dodatku 1 wynika, ze 
w niekthiych przypadkach jest dose duza zgodnosb teorii z obserwacj^. 
Rezultaty obliczeh roznig si^ jednak znaeznie od wynikow opublj'kowa- 
nych przez Bohma w 1929 r. Pewne swiatlo rzuca na powyzsze rozbiez- 
nosci ta okolieznose, ze uj^cie optyki geometrycznej jest sluszne tylko 
dopoty, dopoki cz^stotliwosc robocza nie stanie siQ zbyt bliska granieznej 
czQstotliwosci jonosfery; w tym momencie optyka geometryezna przestaje 
obowi^zywah. 
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Sprawa ta wymaga dodatkowych wyjasnien. Wedlug praw optyki 
geometrycznej przy zmniejszaniu siQ g^stosci jonizacji przejscie od stanu 
odbijania siQ fali jonosferycznej do stand nieodbijania si^ jej odbywa siQ 
skbkowo. Graniczne odbicie zachodzi wtedy, gdy szybkosc grupowa fali 
staje si^ rowna zeru. Pohiewaz jednak iloczyn szybkosci grupowej przez 
fazow^ rowna siQ kwadratowi szybkosci ^wiatla, to w miar^ zblizania siQ 
do czQstotliwosci granicznej (MUF), gdy szybkosc grupowa d^zy do zera, 
szybkosc fazowa d^zy do nieskonczonosci. UwzglQdniajqc jeszcze, ze jono- 
sfera jest niejednorodna, nalezy zauwazyc, ze niejednorodnosci staj^ si^ 
wtedy male w porownaniu z dlugosci^ fali w danym miejscu. W ten spo- 
s6b powstajq wlasnie warunki, w ktorych optyka geometryczna przestaje 
,,obowiqzywac wobec nast^puj^cej ogolnej zasady: metody optyki geome- 
trycznej stosuj^ si^ tylko wtedy, gdy koncentracja elektronowa nie ulega^ 
praktycznie zmianie w granicach. odleglosci porownywalnyeh z dlugosci^ 
fali. ■ ’ 

W rozpatrywanych warunkach musimy siQ zatem uciec dp optyki falo- 
wej. W zwiqzku z tym nizej rozpatrzono mechanizm odbic jonosferycz- 
nych w ujpciu optyki falowej i w oparciu o zdjete doswiadczalnie prze- 
biegi pojedynczych impulsow powracajqcych na ziemip po odbiciu od 
jonosfery. 

3. PRZYPUSZCZALNY MECHANIZM ODBiC OD JONOSFERY 
yv WARUNKACH GRANICZNYCH 

Scisla interpretacja matematyczna zjawisk zachodzqcych w jonosferze 
natrafia w ujpciu optyki falowej na bardzo powazne trudnosci, zwlaszcza 
w poblizu stanow granicznych. • 

Ogolne zagadnienie rozchodzenia sip fal elektromagnetycznych w sro- 
dowisku zjonizowanym znajdujqcym sip w polu magnetycznym ziemi 
nie znalazlo dot^d scislego rozwi^zania. Jest to spowodowane tym, ze: 

1. nie zostaly jeszcze wyprowadzone scisle rownania rozchodzenia sip 
fal w rozwazanych warunkach, 

2. przy rozwi^zywaniu rownah trudno jest sformulowab w sposob 
wlasciwy warunki graniczne. 

Tylko. w niektorych przypadkach szczegolnych rozkladu jonizacji i dla 
niektorych kierunkow specjalnych lub wobec braku pola magnetycznego 
rownania mogp bye sformulowane i moze bye znalezione ich rozwi^zanie. 
Dlatego sprobowano odtworzyb przypuszczalny mechanizm odbic jono- 
sferycznych w tych warunkach w sposob przyblizony, opieraj^c sip na 
pewnych materialach doswiadczalnych. Jako punkt wyjscia rozpatrzono, 
stwierdzone eksperymentalnie, zjawisko rozbicia na szereg ci^gow poje- 
dynezego impulsu, odbitego od warstwy jonosferycznej. 

Przyklad takiego impulsu pokazany jest w Dodatku 3. 
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Mechanizm rozbijania impulsu na ci^gi mozna ewentualnie wytluma- 
czyc kilkoma sposobami, a mianowicie: ' 

1. rozbicie impulsu nast^puje wskutek rozszczepiaj^cego. dzialania 
pola geomagnetycznego w jonosferze, wytwarzaj^cego fal^ zwyczajh^ 
i nadzwyczajn^, 

2. rozbicie impulsu nast^puje wskutek zmiennosci w 'czasie charakte- 
rystyki filtru dolnoprzepustowego, jakim jest jonosfera w poblizu czQsto- 
tliwosci granicznej (o charakterystyce tego filtru b^dzie mowa nizej). 

Rozbity na szereg ciqgow impuls moze bye rozpatrzony za pomoc^ 
calki Fouriera. W przypadku dla przejrzystosci mocno uproszczonym, 
traktuj^c pdszczegolne ci^gi jako odcinki sihusoidy, mozna z latwosci^ 
otrzyrnac na tej drodze rozklad spektralny pojedynezego ci^gu. Przyklad 
takiego rozkladu spektralnego podany jest w Dodatku 4. Jak wynika 

z przedstawionego w Dodatku wykresu dla n=3, wartosci — od- 


powiada wartosc-^^? 

2 


3, CO stanowi okolo 0,64 wartosci maksymalnej- 


Im ci^gi fal s^ dluzsze, tym pasmo staje si^ w^zsze, Dla ci^gow zlozonych 
na przyklad z n=100fal amplituda fall odbitej o cz^stotliwosci rozniqcej 

siQ o — = 0,P/o od czQstotliwosci emitowanej jest prawie rowna amplitu- 
lOn ^ 

dzie fall odbitej o cz^stotliwosci emitowanej; dla cz^stotliwosci emitowa- 
nej 20 MHz stanowi to roznicQ 20 kHz , co jest wartosci^ dose znaezn^. 

Przy stanach bliskich granieznym jonosfera zachowuje si^ selektywnie 
jak filtr dolnoprzepustowy w stosunku do promieni powracajacych na 
ziemiQ. Wynika to z falowego uj^cia zjawiska odbicia, co znajduje swoj 
wyraz w charakterystycznym* przebiegu. wspolczynnika odbicia, ktory 
zobrazowany jest w Dodatku 5. W Dodatku tym, ze wzgl^du na trudnosci 
matematyezne, przedstawiono wspotezynnik odbicia w funkcji roznicy 
czQstotliwosci krytyeznej i biez^cej jedynie dla transmisji prostopadlej 
(czQstotliwosc krytyezna jest szczegolnym przypadkiem czQstotliwosci 
granicznej). Analogicznego przebiegu wspolczynnika odbicia nalezy spo- 
dziewac siQ dla transmisji ukosnej. 

W ujQciu optyki geometryeznej mielibysmy prostok^tny przebieg 
wspolczynnika odbicia R, co nie odpowiada rzeczywistosci, 

Ze wzgl^du na ustawieznie zmieniaj^ce si^ wlasciwosci filtru jonosfe- 
ryeznego i zwi^zane z tym stany przejsciowe, rozklad spektralny impul- 
sow powracajacych na ziemi^ ulega roznym zmianom, zwlaszcza przy kil- 
kakrotnym odbiciu od jonosfery. , 

Jak z tego wynika, szerokie widma wyst^puj^ tylko w warunkach gra- 
nicznych, ktore w praktyce spotyka si^ oczywiscie znaeznie rzadziej niz 
przeci^tne. Dla doswiadczalnego potwierdzenia tego faktu przeprowadzono 
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szereg obserwacji na.radiostacji Instytutu Lacznosci oraz na radiostacji 
odbiorczej Ministerstwa Lacznosci w Grodzisku. 

W Instytucie Lacznosci badania przeprowadzono za pomoc^ nasluchow- 
(z filtrem) i przy uzyciu oscyloskopu. Na.radiostacji odbiorczej w Gro- 
dzisku proby dokonane byly rowniez metod^ nasluchow, a oprocz tego 
przez obserwacji regulowanej automatycznie heterodyny, pod^zajacej za 
czistotliwosci^ odbieranego sygnalu. 

Nalezy zauwazyc, ze opisany mechanizm znacznej zmiany czistotli- 
wosci wyjasnia zwi^zek zachodz^cy miidzy nim a duz^ zmiennosci^ przy- 
chodzenia fal jonosferycznych. 

Reasumuj^c powyzsze rozwazania mozemy stwierdzic, ze przy doko- 
nywaniu pomiarow na podstawie odbioru czistotliwosci wzorcowych od- 
legtych radiostacji krotkofalowych nalezy unikac zblizania siq do stanow 
granicznych jonosfery. 

Instytut L,qczno§ci 
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Dodatek 1 


Dane z literatury odnosnie zjawiska Dopplera 


2r6dlo 

Data 

obserwacji 

Trasa 

i jej diugosc 

. 

Cz^sto- 
tliwosc 
w Mc/s 

— • 10« 
f 

U w a g i • 

Bohm 

1929 r. 

Niemcy- 




[3! 


Buenos 

Aires, 



Z obserwacji fototelegrafii 



ok. 12000 km 

20 

100000 


Both 

1934 r. 

Anglia- 



Dane z obserwacji wzorca 

i Gregory 


Ameryka 



cz^stotliwo§ci, dokladnosc 

[2] 


Polnocna, 



wzorca 



ok. 5050 km 

5 

20 

-108= + 5. , 


Grudzieh 

1946 r. 

Anglia- 

Ameryka 


srednio 

Dokladnosc wzorca w ciqgu 
paru dni 

! 

Godz.9 — 9.15 

Godz. 20.42 

Polnocna, 
ok. 5050 km 

10 

4,7 

' ' 107 

y •10"=± kilka. 


Godz. 9.45 



63 

Zmiany cz^sfotliwoici naj- 


i 



7 

wi^ksze 0 wschodzie i za- 
; chodzie slohca. 

1 Przy zachodzie slohca czq- 
stotliwosc maleje, przy 
wschodzie ro^nie. 


28. 4. 1947— 

Anglia- 


• 

Zmiany cz^stotliwo^ci wi^k- 


—2. 5. 1947 r. 

Ameryka ; 

10 

do 20 

sze dla nizszej czQstotliwosci. 


Godz. 20 — 8 

Polnocna/ 




ok. 5050 km 

15 

10 ‘ 


Griffiths 

Listopad 

Anglia- 



■ 

Dokladnosc cz^stotliwosci 

[4] 

1943 r. 

Ameryka 


2—7 

1 stacji nadawczej 


Godz. 1 . 
Grudzieh 

Polnocna, 
ok. 5050 km 

15 

srednio 3 

+1 

II 

GO 

O 


1945 r. 


15 

3,3 
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Dodatek 2 


gdzie : 


Zjawisko Dopplera przy ruchach zwierciadlanej warstwy jonosfcrycznej 
w uj^ciu optyki geometrycznej 

Ogolny wzor Dopplera 

W zalozenia Vo^ 0 

-f ^ f ' f ' 4, y 

. c + v^-v^ ‘ 

— czQstotliwosc odbierana przez _ c-uco. 

sluchacza' s * / \ ' ^ 

j, — czQstotliwosc wysylana przez / \ 

zrodlo z / \ 

c — pr^dkosc rozchodzenia siQ fal / \ 

w srodowisku / \ 

Vo rzut prQdkosci srodowiska na / ^ 

kierunek zs 

— rzut pr^dkosci sluchacza na 

kierunek zs c ' c-t/cos 2, r. _ 

v„ — rzut pr^dkosci zrodla na kie- c^ucos^^^> c^ucos^2 f 

runek zs 


c c-ucosn, r /^c 

c^ucos^2 f l~^cos^, 


Przypadki szczegolne 


Dla uCc 

f 0 / \l ^ 

( 14 cos/] j 1^1+ — cos 72 


fo ^ 

y^.=y2=90°-a — =1 + 2 — sina 


: n I cos 7i + COS72 



. ^ JU ^ 

a) przy a ^ 0 ; 1 

fo ' ' U Aj u 

b) przy a->90'; >1 + 2 — ; — 

f c f c 



y.x=^ 

y2 = 180°-ajf 


- (cos a — cos a) = 1 


Wg Manczhrskiego [5]. 
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Dodatek 3 


Zjawisko rozbicia na szereg ci^gow pojedynczego impulsu 
odbitego od warstwy 




yiiH 









too 2QO 3CX) 40Q 500 sm^voQ m : 







.200 ..000 40Q. 500 m 


a. Rozbity i rozszerzony sygnal fali 
zwyczajnej; najwi^kszy z ci^gow zaj- 
muje szerokosc ok. 100 km (400 — 300 
km), inne mniejsze ci^gi maj^ szero- 
kosc podobnego rz^du. 


b. Rozbity i rozszerzony sygnal fali 
zwyczajnej, zajmuj^cy szeroko§6 ok. 
220 km (550 -*330 km); ci^g po lewej 
stronie odnosi si^ do fali nadzwy- 
czajnej. 


■ - Mil 



joo 300 '400" . 5aaf^^ , \ 


c. Zespol ci^gow fali zwyczajnej i nad- 
zwyczajnej, zajmuj^cy t^cznie* szero- 
kosc ok. 300 km (500 — 200 km). 


Wedlug zdjqc fotogmficznych, podanych przez Alperta 
w ksiqzce ''Rasprostranienije radiowoln”, 1953 [1]. 
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Zjawisko Dopplera w propagacji 


Rozktad spektralny odcinka sinusoidy 




S — modut g^stosci spektralnej, * 

CO — pulsacja dowolnej skladowej w pasmie, 
coq— pulsacja czQstotliwosci emitowanej, 
n — ilosc fal sinusoidy. 

0 najblizej max. dla i-— = — 

2 I a)o 1 n 


Dodatek 4 


Dodatek 5 


Obliczony wspolczynnik odbicia od warstwy F w uj^ciu optyki falowej 




I 




\Ar^{ro-f)c/s 


-200 -wo 0 too 200 300 ^00 500 600 700 800 

R — wspolczynnik odbicia,’ 

V — ilosc zderzen na sekundQ elektronow z molekulami, 
fo — cz^stotliwosc krytyczna, 
f — czQstotliwosc fali. 

Obliczenia odpowiadaj^ zalozeniom : 

a) paraboliczny rozklad gQstosci jonizacji w warstwie, 

b) polowiczna grubosc warstwy — 120 km, 

c) fo - 10 Mc/s. 

Wg Ginzburga, '’Rasprostranienije radiowoln”, 1953 [IJ, [71. 
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C. BOPOBCKH, C. 5iCHHbCKH, C. MflHHflFCKH 

HBJIEHME ^OnnjIEPA B MOHOC<^>EPHOM PACIIPOCTPAHEHMM 

PA^MOBOJIH 

P e 3 K) M e 

Ii;ejibio HacTOflmeM paGoTbi HejiHCTCH MS'bHCHeHMe ECTpenaeMbix ,b jiMTeparype 
B033peHHM Ha TCiviy npOMCxoAHmero, b cbhsm c HecTaSmibHOCTio MOHOcct)epbi, hbjic- 
HPia ^onnjiepa m ero 3HaMeHna b kopotkobojihobom paAMOCBaaw. UjpiTMpoBaHHbie 
pasHbiMM aBTopaMM pesy.nbTaTbi wsMepenpipi KOJie6aHMa aacTOTbi, Bbi3biBaeMoro 3tmm 
aBjieHMeM (npwjiOHceHMe 1), OTJiHaaiOTca APyr or T^pyra na HecKOJibKo fiopa^KOB. 

HeKOTOpbie aBTopbi Bw^aT npMHMHbi B03HHKaHMa CMJibHoro HBJieHMH ^onnjiepa 
B GOJIbmOM MSMCHHHBOCTM yrJIE npMXO^a MpHOC(|)epHOfi BOaHbl. 

;];oKJiaA co^[^epatnT Kax MaTeMaTMHecKMM anajiMB npohJicMbi (npHJioxceHMe 2), rax 
II pesyabTaTbi wsMepeHMif m Ha6jiK);iieHMM, MareMaTiiMecKaa (^JopMyjia KoaeSaHwii 
nacTOTbi BcaeACTBMe hoabiixchoctm MOHOC<J)epbi (npnaoxceHMe 2), BbiBeAennaa Ha 
ocHOBax reoMeTpMHecKOM onTMKM, Aaer, Kax. aro cjieAyer ns BbiHMcaeHMM, Heana- 
HMTeJIbHbie M3MeHMH HaCTOTbl. 3 tO. O^HaKO, npaBMJIbHO TOJIbKO AO MOMCHTa, KOrAa 
paSonaa nacroTa ne craHCT cjimuikom Gjiiiskoh k npeAOJibHOM HacTOxe woHOccJ^epbi 
MIIM. B 3TOT MOMCHT npeKpamacTCH npMMeHeHMe reoMeTpMHecKOM onTMKM: npM 
cocToaHMHx" 6.JIM3KI1X K npcAeabHOMy MsSiipaTejibHoe noBeAenwe woHoccJ^epbi hoaoSho 
cJjHJibTpy hmxchmx aacTOT no OTHOmeHMio k B02Bpamaiou4MMca Ha aeMJiio JiyaaM, 
npMHeM CBOMCTBa 3TOrO (|)MJlbTpa nOCTOaHHO HSMCHHIOTCa. B CJieACTBMe B03HM- 
KaioiHHx TaKnM oGpaaoM nepexoAHbix' coctohhmm BOBBpamaioinMMca na aeMJiio, 
OTAejibHbiii MMnyabc noABeprae'Pca pacmenjiCHmo na nocjieAOBaTeabHbie KOpOTKwe 
rpynnbi (npMaojKeHiie 3). Ecjin aareivi 3 tw rpynnbi paccMarpiiBaTb npw noMou^n 
MHTerpajia Oypbe (npiiaoxceHMe 4), to noayHaeTca aobojibho inHpOKiiM cnexTp nacTOT. 
npM6jm3HTeJIbHO COOTBeTCTSyiOIUMM A^BIiblM HCKOTOpblX aBTOpOB. 

PacnpeABJieHiie aacTOT b cnexTpe BosBpamaiomeroca na seMJiio cHrnajia Moxcex 
6biTb paaawHHO (npwji. 4, 5). Ms SToro cjieAyex, hto ninpoKiie cnexTpbi BbicxynaioT Toab- 
Ko B npcAeabHbix ycjiOBnax, KOTopbie Ha npaKTMKe BCTpenaiOTca oneBMAHO ropasAo 
pex^e HCM cpeAHwe. ^jih onbiTHoro noATsepacAOHHH 3Toro (JiaxTa 6bia McnojineH paA 
Ha6jiioAOHiiM paAMOCTaHi^Mew MHCTHTyxa Cbhsm,’ a xaKace npiieMHOM cxanqHeii Mh- 
HHCxepcxBa CBa3n b PpoABiicKe. CacAyex saMexHXb, hxo onwcaHHbiii MexaHMSM cMab- 
Horo aBaeHMH T^onnaepa xopomo MSTbacHaex cBa3b>Me5KAy othm HBJieHMeM w ana-' 
HMxejibHOM M3MeHHMB0CTbK> yraa npifxoAa iiOHOC(J)epHbix bojih. 

BbixeKaiouj^MH m 3 Hacxoari^eM pa6oTbi npaKXMHecKMH bbiboa tox, hxo npM ricnoji- 

KeHMH M3MepeHMM Ha OCHOBaHMM HOpMaJIfaHbIX HHCTOT OTAaJieHHbIX KOpOTKOBOaHOBbIX 

cxanquM — caeAyex ii36eraTb opnGaiixceHMa k npCAeabHbiM cocxoaHHHivi M0H0ccJ)epbi. 

S. BOROWSKI S. JASmSKI, S.‘ MANCZARSKI 

DOPPLER EFFECT IN IONOSPHERIC PROPAGATION 
Summary 

The object of this paper is to explain the difference in views found in the lite- 
rature on the importance of the Doppler effect in short-wave communications, this 
effect arising from the variations of the ionosphere. Measurements of frequency vara- 
tions due to Doppler effect, given by various authors, vary widely, even to the extent 
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of several orders (Appendix 1). Some authors associate the appearance of a strong 
Doppler effect with the great variability of the angle of incidence of the waves on 
the ionosphere. 

This paper contains a mathematical analysis of the general problem (Appendix 2) 
as well as an analysis of the measurements and observations of the effect. A formula 
expressing the variation of frequency caused by the motion of the ionosphere 
(Appendix 2) has been derived on the basis of geometrical optics. Calculations made 
using this formula give a small change of frequency. The formula can be used as 
long as the operating frequency is not too close to the ionospheric MUF fre- 
quency. In this case geometrical optics is no longer applicable. Near this limit the 
ionosphere behaves like a selective low-pass filter for the waves returning to the 
earth, whereby the parameters of this filter are constantly changing. As a result of 
the changes arising in this manner, a single impulse returning to the earth is broken 
up into a series of short wave trains (Appendix 3). If these wave trains are studied 
with the aid of the Fourier integral (Appendix 4), we obtain a wide frequency spec- 
trum, which more or less is in agreement with the data obtained by some .authors. 
The spectral frequency distribution for the signal returning to the earth can vary 
widely (Appendix 4 and 5). . 

Analysis shows that the wide spectra appear only in the limiting conditions, 
which certainly occur considerably less often than the average ones. To confirm 
this fact experimentally, a series of observations, were made at Warsaw and Gro- 
dzisk (Poland). It should be noted that the treatment presented by the authors pro- 
vides a good explication of the correlation between the Doppler effect and the great 
variability of^ the angle of arrival of ionospheric waves. 

As a practical conclusion results, it may be noted, that in carrying out mesur- 
ements based on the. reception of standard frequencies from distant short-wave ra- 
diostations, an approach to the limiting conditions should be avoided. 


9 ‘ ■ 
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621.3.094.3:621.396,82 

R. KULIKOWSKI, A. RYBARSKI 

Kryterium znieksztalcen kwadratowych 
przy ograniczonej mocy zaklocen 

R^kopis dostarczono 31. 10. 1955 


Wyznaczono optymalne charakterystyki ukladu liniowego zapewniaj^ce 
przeniesienie sygnalu o zadanym ksztalcie z minimalnymi znieksztalceniami 
kwadratowymi przy zalozeniu, ze srednia moc zaklocen jest ograniczona. 
Wyznaczono sygnaly pozyteczne i zaklocaj^ce, ktore najbardziej nieko- 
rzystne z punktu widzenia pracy ukladu. Rozpatrzono rowniez sygnaly modu- 
lowane i charakterystyki ukladow selekty wnych. 

\ , 

1. WSTIJP 

W teorii odbioru sygnalow elektrycznych duzo uwagi poswi^ca siQ 
tzw. filtrom liniowym, ktorych zadaniem jest wydzielanie sygnalu po- 
' zytecznego sposrod zaklocen. Dla okreslenia i porownania zdolnosci fil- 
truj^cych filtrow odbiorczych liniowych stosuje si^ roznorodne kryteria. 
Wybor kryterium uzaleznia siQ zwykle od konkretnych warunkdw eks- 
ploatacyjnych, w jakich pracuje urz^dzenie odbiorcze. 

Jednym z najbardziej rozpowszechnionych kryteriow jest stosunek 
maksymalnej mocy sygnalu do mocy zaklocen (lub stosunek amplitudy 
sygnalu do- skutecznej-amplitudy . zaklocen) ,[!]. 

Oznaczaj^c przez S(ift>) widmo sygnalu, a przez Zi(ift>) widmo zaklocen, 
kryterium powyzsze mozendy wyrazic jako 

+ 00 

2 1 J iS(ico)K(ift>) | ^ 

= , ( 1 ) 

• / 1 Z(i(o)K(ico) I ^dco 

— oo 

gdzie K(ico) jest funkcj^ przenoszenia ukladu, T zas symbolizuje moment 
v^czasu, dla.ktorego amplituda sygnalu osi^ga wartpsc-maksymaln^. 

^ Korzystaj^c z nierownosci Cauchy-Buniakowskiego mozna wykazac, 
ze stosunek powy^zy osi^ga maksymaln^ .wartosc dla filtru .o funkcji 
przenoszenia. 
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K(m) = 




gdzie oznacza funkcj^ zespolono-sprz^zona z S(ico). 

Kryterium powyzsze okazuje siq szczegolnie wygodne w przypadku, 
gdy sygnal *S sklada siq z szeregu identycznych impulsow przedzielonych . 
znacznymi (w porownaniu z szerokosci^ impulsu) odcinkami czasu (jak 
to ma miejsce na przyldad w radiolokacji). W przypadku tym mozna rze- 
czywisty sygnal zast^pic' pojedynczym impulsem, czyli inaczej mowi^c, 
okres powtarzania impulsow rozci^ga si^ do nieskonczonosci. 

Nalezy podkreslic, ze kryterium stosunku sygnahi do zaklocen nadaje 
siQ szczegolnie do ukladow, w ktorych zalezy glownie na wykryciu sygna- 
!6w (b^dz momentow ich pojawienia si^), natomiast nie zalezy na doklad- 
nym odtworzeniu ksztaltow sygnalu. W tych urz^'dzeniach, gdzie wazne 
jest mozliwie dokladne odtworzenie ksztaltu sygnalu, szerokie rozpo- 
wszechnienie uzyskalo kryterium dparte na tzw. bl^dzie sredniokwadra- 
towym [5]: 


gdzie 


7 = lim — r [s{t)--u(t)fdt , 
T->oo 2T J 

~T 


CX) 


t) k(r) dr , 


e(t) — s(t)+z(t) — jest sygnalem na wejsciu ukladu, 
s(t) — sygnalem pozytecznym, ; ’ 
z(t) — sygnaleni . zaklocenia, 

k(t) — funkcj^.przejsciow^, ktora jest reakcj^' ukladu na 
sygnal impulsowy 5(t). 

Jesli zalozymy, ze sygnal i zaklocenia opisuj^ siQ funkcjami stocha- 
stycznymi stacjonarnymi (w szerokim sensie), to okazuje si^, ze wyzna- 
czenie opfymalnej funkcji przejsciowej k(t) sprowadza siQ do rozwi^zania 
rownania calkowego typu Wienera-Hopfa 


gdzie 


Res(t)— I k(r) Re(t—r) dr^^O , 


Res=lim — ( s(t— t) [s(T)+2:(r)] dr , 
T^oc2Tj 

— T 
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Kryterium znieksztalcen kwadratowych 


5 =lim — r 

T-^ 2TJ 


e(t— t) e(r) dr 


sq tak zwanymi funkcjami korelacyjnymi sygnalu s{t) i z{t). 

Wynika st^d, ze w zagadnieniu filtracji, opartym na kryterium (3), 
wystarczy znac funkcje korelacyjne sygnalu i zaklocen. Z geometrycz- 
nego punktu widzenia zagadnienie sprowadza si^ do przyblizenia funkcji 
s(t) i u{t) w Hilbertowskiej ^przestrzeni funkcyjnej, ktorej wlasnosci okre- 
sla iioczyn skalarny wektorow s{t) i u{t), iloczyn ten pokrywa siQ nato- 
miast z funkcje korelacyjn^. 

W kryteriach (1) i.(3) operuje siQ b^dz (Stosunkiem mocy, b^dz znie- 
ksztalceniami sygnalu. Kryteria te nie uwzgl^dniaj^ wi^c wszystkich 
czynnikow interesuj^cych w praktyce. Na przyklad nie uwzglQdniaj^ one 
wzajemnego stosunku mi^dzy moc^ (lub energi^) sygnalu i zaklocen na 
wyjsciu ukladu odbiorczego, a znieksztalceniami sygnalu.' 

Dlatego tez wydaje si^ celowe zbadanie kryteriow opartych z jednej . 
strony na znieksztalceniu sygnalu, z drugiej — operujqcych moc^ (lub 
energi^) sygnalu oraz zaklocen. 

W niniejszej pracy ograniczamy si^ do rozpatrzenia kryteriow, w kto- 
rych dla oceny znieksztalcen wykorzystuje si^ bl^d kwadratowy przy 
ograniczonej mocy (lub energii) zaklocen na wyjsciu ukladu. 


2. KRYTERIUM ZNIEKSZTAECElSr KWADRATOWYCH 
PRZY OGRANICZONEJ MOCY ZAKLCCEN 

Gdy sygnal e(t) dziala na uklad liniowy o funkcji przejscia h{t) wyra- 
zaj^cej reakcjQ ukladu na skok jednostkowy 


lOdlaKO, 

wowczas na wyjsciu ukladu ukaze siQ sygnal 


u{t) = — I e(t — r)h(r)dT . 
dt .1 


W rozmaitych zagadnieniach wymaga si^, aby sygnal wyjsciowy z mi- 
nimalnym znieksztalceniem odtwarzal przebieg x{t). Na przyklad przy 
projektowaniu ukladow przekazuj^cych i wzmacniaj^cych mamy zalez- 
nosc •• 

x(t) = const ■ e(t— to) 
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gdzie to jest zadanym czasem opoznienia sygnalu w ukladzie. Jako kry- 
terium sluz^ce do oceny znieksztalcen rozwazymy wyrazenie 

OO ' 

' 1= j‘[x{t)-u(t)fdt. (6) 

0 

Oprocz sygnalow pozytecznych na uklad dzialaj^ czQsto sygnaly zaklo- 
caj^ce, ktorych efekt dzialania powinien bye dostateeznie wygaszony, 

Efekt dzialania zaklocen mozna oceniac roznie. Tu ograniezymy siQ 
do rozpatrzenia sredniej mocy szumow oraz energii zaklocen na wyjsciu 
ukladu. 

Jesli na przyklad g^stosc .widmowa mocy szumow wejsciowych moze 
bye uwazana w przedziale cz^stotliwosci obj^tych charakterystyk%.filtru 
za wartosc stal^ //, to srednia moc szumow wyjsciowych wyniesie 

+ 00 

— r I K(i(o) \^d(o . 

' j • - 

— OO 

Mozna wykazac, ze dla wszystkich zaklocen z ograniezon^ powierzch- 
rii^ /< energia zaklocen wyjsciowych nie przekraeza wartosci okreslonej 
powyzszym wzorem. Celowosc operowania powierzchni^ zaklocen zwi^- 
zana jest z tym, ze przy dzialaniu na uklad liniowy krotkich impulsow 
z(t), reakeja jego prawie zupelnie nie zalezy od ksztaltu sygnalu, lecz je- 
dynie od jego powierzchni. 

Wyznaczymy teraz sposrod wszystkich sygnalow zaklocaj^cych o ogra- 
niezonej powierzchni 

- +00 

J| 2(t) I dt = /i, (7) 

— OO 

takie, ktorym odpowiada maksymalna energia zakldcen na wyjsciu ukladu 


Tutaj 


2icJ 


Z{ico)= J 2 (t)e“*‘“‘dt , K(ico)= J fc(t)e“‘“‘dt . 


Wykazemy, ze sq to sygnaly o stalej g^stosci widmowej. Dowod opiera siQ 


na nierownosci 


I Z(icoj 1 < J i z(t) I • I e-‘"* ! dt= f | z(t) ldt=ju, 
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sk^d dla energii zakiocen otrzymujemy oszacowanie 

^ +00 -foo . 

E= — r |ZK|"dft)<^ r I K(i<M) I const. 

2x j ' ' 2 jc j 

— 00 ~O0 J 

Wobec tego najwi^ksza wartosc energii E jest osi^gana dla sygnalow, przy 
ktorych obie nierownosci przechodz^ w rownosc, skqd otrzymujemy 

\Z{i(a)\= [j,. 

Sygnalem takiego typu je§t na przyklad kombinacja liniowa impulsow 

k ' 

poprzesuwanych w f azie. 

Dla okreslonych powyzej najniekorzystniejszych zakiocen energia na 
wyjsciu wyrazi siq wzorem 

+00 

E= f \K{i(o)\Ha,. (9) 

—00 

Obecnie sformulujemy scisle rozpatrywany problem. 

Niech b^dzi'e zadany sygnal wejsciowy e(t) . Nalezy znalezc uklad 
odtwarzaj^cy z minimalnym odksztalceniem (6) zadany przebieg.x(t) , przy 
' czym energia zakiocen na wyjsciu nie moze przekraczac zadanej wartosci 
E. Jezeli zauwazymy, ze z pomocq zwi^zku h'{t)=k{t) wi^zqcego wprowa- 
dzone funkcje przejscia, na mocy rownosci Plancherera mozemy zapisac 
energi^ zakiocen wzorem . 

oo 

Jh’^(t)dt, (10) 

0 

to zadanie nasze, z punktu widzenia matematycznego, jest zadaniem z ra- 
chunku wariacyjnego na minimum funkcjonalu (6) przy izoperymetrycz- 
nym warunku'ubocziiym (10). 

Przy podobnym podejsciu do zagadnienia mozna wyznaczyc opty- 
malnq charakterystyk^ ukladu poddanego dzialaniu sygnalu niesto- 

chastycznego o zadanym ksztalcie. Jesli zas sygnal dzialajqcy na uklad 
sklada siQ ze skladowej stochastycznej stacjonarnej e*, niestochastycznej Cn, 
i zakiocen stochastycznych stacjonamych Zg, to mozna wykazac [6], ze 
srednia reakcja ukladu nie zalezy od skladowych stochastycznych, lecz 
okresla si^ skladow^ niestochastycznq e„ . Mozna wi^c znalezc minimum 
bl^du kwadratowego (6) I dla przebiegow niestochastycznych przy zalo- 
zonej wielko^ci bl^du ^redniokwadratowego (3) I = const (w ktorynx 
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uwzgl^dnia si^ skiadowe stacjoname) oraz przy zalozanej sredniej mocy 
zaklocen E=const (lub tez minimum I przy I=const, E=const). 

Nieco irmq metodQ podano w pracy [6], gdzie sktadow^ uiestocha- 
styczn^ rozwija si^ w szereg w przedziale [0, T] 

en(t—r) = en{t)- e'Jt)r + — 

i otrzymuje si^ odpowiednio 

n ' 

)c=l 

• gdzie 

T 

' fiii= j r’^k(r)dr , fc=0,l,2, n. 

0 


Przyrownujqc nast^pnie u„(t) do zadanej funkcji czasu otrzymuje si^ 
warunki na f^k. Wyznaczenie optymalnej funkcji przejsciowej k(t) spro- 
wadza si^ tu do znalezienia minimakiej wartosci funkcjonalu I przy wa- 
runkach dodatkowych nalozonych na /ik • Nie trudno zauwazyc, ze metoda 
powyzsza komplikuje siQ w przypadku bardziej zlozonych sygnalow (na 
przyklad sygnalow, ktore nie mog^ bye przedstawione wielomianem alge- 
braieznym dostateeznie niskiego rzQdu). 

Poniewaz zagadnienia zwi^zane z okresleniem optymalnych charakte- 
rystyk ukladu przy sygnalach stochastycznych i stacjonarnych sa wyezer- 
puj^co omowione w literaturze, zajmiemy siQ ponizej, glownie wyznacze- 
niem charakterystyk ukladow poddanych dzialaniu sygnalow niestaejo- 
narnych o zadanym ksztalcie. 



2.1. ROWNANIE CHARAKTERYSTYK OPTYMALNYCH 

Zgodnie z metod^ . mnoznikow Lagrange’a rownanie na funkcji h(t) 
znajdziemy przyrownuj^c do zera wariaej^ funkcjonalu pomocniczego 

W=r + AE', (11) 

w ktorym A oznaeza mnoznik Lagrange’a. 

Jesli r](t) b^dzie dowolnq funkcji znikajqc^ poza skonezonym przedzia- 
lem i rown^ zeru dla t=0, wtedy warunek ekstremum funkcjonalu (11) 
mozemy napisac jako 

aW=l— W[h+e»?]} =0. (12) 

I ds J a=0 

W dalszym ci^gu obliczymy dW . Mamy 

2 ds 2 de 2 de 
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Ze wzorow (6) i (10) po przeksztalceniach otrzymamy 

oo oo 

~=“ diri(t) |'d,e(T-t)[a:(T) -u(T)], 


dl 

2 


2 d£ J 


dt , 


a- nastQpnie 


' oo oo 

dw= — J" ?7(t)dt| Ah"+ J* d^ 


?7(t)dt [/h"+ j d^e(T— t) [ar(T) — u(t)]| . 


Obecnie z warunku ekstremum (12), wykorzystuj^c dowolnosc funkcji r/(t) 
otrzyrpujemy rownanie 


Ah"(t)-h J d^e(T— t)[x(T) — u(t)] = 0. 


(13) 


Rownanie to musi spelnia6 optymalna charakterystyka-h(t) . Przepiszemy 
je nieco inaczej, przyjmuj^c dla skrotu 

OO ;■ oo 

b(t)= Jde((j)e((T+t); c(t)= j'de(a)x(a + t). (14): 

. 0 ' 0 " 

Wtedy 

' I ‘ oc tx) » 

J d,e(r - 1) u(t) = J dxh'{r) b(t - r) 

0 0 

i rownanie (13)’ mozemy przepisac w postaci 

oo 

~^h"{t)+ jb(t—r)hXr)dT=c{t). (15), 

- ® , 

Obecnie iidow^odnimy, ze dla rozwi^zanie rownania (15) realizuje 

minimum funkcjonalu (11). W istocie, obliczaj^c drug^ wariacj^ d^W 
ze wzoru 




znajdujemy 


Ide® 


W[h + srj] 


c=0 


oo . 00 
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u— — e{t — r)[h(r) + er]{r)]dr. 
dtj 
0 


Archi Elek'trot. 


Poniewaz A<0, wi^c 5W>0, jak wiadomo,- jest wystarczaj^cym wa- 
runkiem minimum. - , . 

2.2. JEDNOZNACZNOSC ROZWIAZANIA ROWNANIA (15) 

Gdy k(t) — h'(t), wtedy zamiast (15) nlozemy napisac 



• Afc'(t) + j b (t — t) k(r) dr = c(t) 


Rownania (15) i (16) mozna sprowadzic do rownan czysto calkowych, . wy- 
korzystujgc rownosci 


k'(t)= J d'{t—r)k(r) dr—k(0) d{t ) , 
0 

oo 

b"(t)= J d"(t—ryh(r) dr — k(0) d(t) . 


Oznaczaj^c 


I,(t) = -Ad'(t) + b(t), I,(t)=-^d"(t) + b'(t), 

zamiast (15) i (16) otrzymamy rownania 

OC 

J Ii(t — T)k(r) dT = c(t)~^k(0) d(t) J (15a) 

0 

oo 

Jl 2 (t-T)h(f)dT+[b(+ 0 )-b(- 0 )]b(t) = c(t)-Ak(0)(5(t). . (16a) 

0 

Rownania te sa odpowiednio typu Focka i Wienera-Hopfa [3], [4]. W pra- 
cach tych autorow sformulo\yane sa zalozenia, przy ktorych rownania takie 
przy zadanej prawej strohia maja jedno i tylko jedno rozwiazanie. 

Prawe strony rownan (15) i (16) zawieraja niOwiadoma wielkosc ?c(0) , 
a wiQC twierdzen tych nie mozemy stosowac bezposrednip. Trzeba jednak 
zauwazyc, ze z uwagi na warunek izoperymetryczny (10) interesuja nas 

oo ^ 

jedynie rozwiazania, dla ktorych J k^dt < + oo . Otoz mozna wykaza6, ze 

0 

oo 

zbieznosc tej calki (i calki J kdt) jednoznacznie wyznacza .k(0) .• Wobec 

.0 . 

tego mozemy stwierdzic, ze nasze rownania iilaj^ co hajwyzej jedno roz- 
wi^zaniQ w interesuj^cej nas klasie funkcii (mozna wykazac, ze wynika 
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to rowniez z rownania (12), bez powolywania sie na teori^ Wienera- 
Hopfa). 


2.3. UKtADY PRZEKAZUJACE 

Zajmiemy si^ teraz przypadkiem, gdy a:(t)=e(t-ro) , to jest gdy 
sygnal, ktory, z minimaniym znieksztakeniem ma bye odtworzony na 
wyjsciu, jest sygnaiem pierwotnym, opoznionym o t„ . Jak to wynika ze 
.wzorii (14), mamy wtedy c(t)=b(t-to) , a rownanie (16) przyjmie postac 

oc 

— Afc'(t)+ f b(t—T)k(T)'dr=b(t—to) . (17) 

0 

oo 

Rozwi^zanie tego rownania, spelniajace warunek Jfc^dt< + oo, jest funkcj^ 

0 

parametru to . Zamiast (17) mozemy wiQc napisac 

CXI 

-Ji^k(t,to)+ J b(t-T)k(r,to)dr=b(t-to). (18) 

0 

Niech /Co(t) = k (t , 0) . Zatem fco(0 j^st optymaln^ charakterystykq, odpo- 
wiadajac^ czasowi opoznienia 0. Wykazemy, ze charakterystyki k(t,to) 
i ,ko(t) zwi^zane s^ rowhaniem 

0 0 

--fc(t,to)+ — fc(t,to) + k(t,0)-?c(to,0) = 0 (19) 

ot 0to 

Dla dowodu napiszmy rownanie (18) raz dla to = ti, drugi raz dla to = t 2 ; 
pierwsze z tych rownan pomnozmy przez k(t , ta) , drugie zas przez k{t , tj) ’ 
a nastQpnie dodajmy i scalkujmy po t od zera do nieskonczonosci. Otrzy- 
mamy 


cx oo 

;.k(0 , tj)fc(0 , ta) + yJ' J dtdT|b(t - t) + b(T - 1)| |fc(T, t j)fc(t , t^) + k(t , t2)k(T , t, )| = 


oo cx 

.= J b(t-ti) k(t,t,)dt+ fb(t-t,)k(t,t,) dt . 
0 0 


Nast^pnie wykorzystamy wynikajqc^ z (14) rownosc i 

'■b(t) + b(— t) = e^( + oo)=e^oc 
oraz fakt, ze, jak mozna udowodnic, mamy 



Jk(t,to)dt=l dla 

0 , . ^ 

1 Przy oczywistym zalozeniu e(t)=0, dla ;t<d. 
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Wtedy ptrzymamy 

OO OO 

;k(0,t,)fc(0,t,) + e^ J k{t,t,)dt f k(t,.t 2 )dt--= 

0 0 

oc • OO ^ » OO OO 

= eifk (t.ta) dt-fb (t,-t)k(t,t^) dt-he^fk (t,t^^dt-Jb (tg-t) Jc dt. 

0 0 0 0 

Calki z prawej strony mozna wyrugowac za pomoc^ (18), przy czym skrocq 
si^ wyrazy zawieraj^ce e^. Dziel^c przez A otrzymujemy 

— k (t, , g +--k (t2 , h) + fc (0 , ti) k (0 , tg) = 0 . (20) 

ati dti 

Mozna udowodnic (przy dodatkowych zalozeniach o zachowaniu sip k(t) 
w nieskonczonosci), ze k(t,to) jest symetrycznp funkcjp swoich argu- 
mentow 

k(t,to) = k(to,t) . (21) 

Ze zwipzkow (20) i (21) wynika (19). 

Rownanie (19) pozwala wyznaczyc k(t,to), o ile znamy ko(t). Aby to 
wykazac, poddajemy (19) transformacji Laplace’a wzgledem t i to , ozna- 
czajpc zmienrip operatorowp odpowiadajpcp t przez p, a zmiennp opera- 
torowp odpowiadaj^cp to przez q.. Oznaczajqc dalej 

|k(t , to)} - K(p, q) ; J?t{ko(t) 1 = Ko(p) I A {ko(to)) = Ko(q) 

w ten sposob z (19) otrzymamy rownosc 

, Ko(p) + Ko(q)-Ko(p)-Ko(q) 

K(p,q)= 7 • 

p + q 

Przez wzipcie transformatory odwrotnej wzglpdem q, dostajemy stqd, jesli 
oznaczymy K(p,to) = j2{k(t,to)} , rownanie 

K(p,to) = Ko(p)e-P‘» + [l-Ko(p)]e-P‘“/ko(T)eP'dT. (22) 

0 

Rownanie (22) pozwala znalezc k(t,to) jesli znamy fc(t, 0) = ko(t). 


2.4. CHARAKTERYSTYKI BEZ OP02NIENIA 

Zajmiemy sip rownaniem (18) dla zerowego czasu opoznienia to=0. 
Wtedy 

OO 

-.lk;,(t)+ / b(t-T)ko(T)dT=b(t). (23) 

0 

Jak juz pokazano, rownanie to mozna sprowadzic do rownania typu 
Wienera-Hopfa, dla ktorego podane sq ogolne wzory na rozwiqzanie. Tutaj 


i 


4 
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pokazemy, ze gdy transformata Laplace’a funkcji b(t) , B(p) jest funkcjq 
wymiern^, to istnieje rozwiazanie o postaci 

V A -f 


ko(t) = 




czyli Ko.(p) 


p-ai 


ou 

Powinnismy miec / fcjdt<oo. Wynika st^d, ze liczby Uf musz^ miec 
0 

ujemne cz^sci rzeczywiste. Wyznaczymy liczby ai iAi, Podstawiaj^c (24) 
do rownania (23) i transfofmuj^c otrzyman^ rownosc wedlug Laplace’a, 
po uporz^dkowaniu otrzymamy rownanie 

'S^—-l-^ai+B(p)+Vi]=B(p)', (25) 


oo 

,/m- 


a) e^i^da = 


-B(-aO. 


Mnozqc rownanie (25) przez p — a^ i zakladaj^c p= otrzymujemy 

-/lafc + B(ofc) + Ffc=0 


-Aafc + B(a;e)- — -B(-afc) = 0 . 

A wiQC liczby ai pierwiastkami o ujemnej cz^sci rzeczywistej rownania 
chardkterystycznego 

e' 

^p+^+B(-p)-B(p) = 0. (26) 

Aby teraz wyznaczyc liczby Ai, porownujemy residua lewej i prawej strony 
(25). Jesli przez oznaczymy bieguny funkcji B(p) (a mozna wj^kazac, ze 
jest ich akurat tyle, ile ujemnych pierwiastkow rownania 26), to dosta- 
niemy 

V 

jezeli biegun Pk jest pojedynczy oraz ' 

2’, 7^=0, , . ^ (27b) 

dla bieguna o krotnosci Tifc> 1 . 

Powyzsze rownania ze wzgl^du na Ai mozna rozwi^zac w postaci 
zarnkni^tej poslugujqc si^ pomocniczymi funkcjanoi 

F(p)=n{p-ai); G(p)=/7(p-|fffc)"*. (28) 

i k 
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Mamy mianowicie 


dla Ko(p) zas 


F'iai) ’ 


Ko(p) = l- 


Maj^c Ko(p) znajdujemy K(p, to) ze wzoru (22). 

Na koniec mozna wykazac, ze przy ai i Ai wyznaczonych z (26) i (29). 
rownanie (25) jest tozsamosciowo speinione, a tym samym wzory (24) 
daj^ rzeczywiscie rozwi^zanie rownania (23). 


2.5. PRZYK£.ADY ’ 

% 

Dla przykladu znajdziemy optymalne charakterystyki dla sygnalu 
\ . 

^(t)=e“‘, x{t) = e(t—to), 

Ze wzoru (14) otrzymujemy 

. B(P)= 

2 2 l+p 

Rownanie charakterystyczne (26) ma postac 

1 — p^ 

i jeden pierwiastek ujemny a= — ^ obliczamy pomocnicze 

funkcje F i G jako 


W koncu z (29) otrzymujemy 


i z (24) 


=p+|/l + y: G = p + 1. 

ijemy 


ko (t)=Ae"‘ = ( + ^ - l) 


Ko(p) = 




P-“ P + 
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7\Vz6r ten podaje optymalne charakterystyki zerowe. Nast^pnie z (22) znaj- 
dujemy 


K(p,to) = 


■A f'at 


p2— a® 


e“‘“(p + l) + e 


A e“‘“ I ch at H — ^ sh atj dla , 


A ( ch atj. H sh atg ) dla t > to 


/c(t;to) = 


to charakterystyki optymalne dla czasu opoznienia to . 

Jako drugi przyklad wyznaczymy charakterystyki optymalne dla 
sygnalu e(t) =l(t), x(t) = 1 (t-to) . 

Ze wzoru (14) mamy 

b(t) = l(t), 

zatem 

£{b(t)}=B(p) = -. 

P 

Rownanie charakterystyczne (26) ma ksztatt nastepuj^cy 


oraz pierwiastek ' 


;p__=0, 

p 




Zatem fuhkcje pomocnicze (28) przyjmuj^ postac 

F=p—a, G—p, 


■ k V i 


Korzystaj^c ze wzorow (30) i (22) otrzymujemy 


Ko(p) = 


p + A 


p^- A^ e-^t»- A e-pt«| 


Wobec tego 


k{t, to) = 


|Ae--*‘»chAt, t<to 


Ue-'*‘ch Ato, t>to. 

Charakterystyka o identycznej postaci zostala' wyznaczona w pracy [2]. 
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Charakterystyka odpowiadaj^ca sygnalowi jednostkowemu jest bardzo 
interesuj^ca z praktycznego punktu widzenia, gdyz — jak to zaraz wyka- 
zemy — skok jednostkowy jest sygnalem najbardziej niekorzystnym 
z punktu widzenia zhieksztaicen sygnalu (w metryce przestrzeni C) z po- 
miQdzy wszystkich sygnalow z ograniczon^ amplitude. 

W tym celu zalozmy, ze sygnal e{t) dziaia na uklad liniowy, ktorego 
reakcja na impuls jednostkowy jest k{t) . Sygnat na wyj^ciu uktadu wy- 


u(t)= j e(t — r)k (t) dr . (32X 

0 

W przypadku, gdy uklad nie wprowadza znieksztaicen 

= (33) 

to znaczy, ze odpowiedz ukladu na sygnai d (t) rozni siq od tego sygnalu 
tylko opoznieniem to i tlumieniem fco . 

Podstawiaj^c (33) do (32) otrzymujemy 

CXI • 

u(t)= Je(t—T)koS (T—t„)dr = koe (t—to) . (34) 

0 

W przypadku ukladu, ktory ' wprowadza znieksztalcenia mamy 
u(t) ^koe(t~to) , Znieksztalcenia te mozna wyrazic za pomoc^' wzoru 

Zl“max|u(t)“koe(t“^to)| . (35) 

t 

Podstawiaj^c wzory (32), (34) do (35) mozemy napisac 

oo oo 

•/l=max j J e(t — r) k (r)dr — j fc^e (t — t) <5 (t — to) dr | = 


max I J e(t — t) [k (r)-~ka<5 (t — to)] dr | . 


Sygnaly dzialajqce na uklad posiadaj^ zwykle ograniczon^ amplitude i dla- 
tego mozemy zalozyc 

|e(t)|<E^. . _ 

Wyznaczmy teraz sposrod sygnalow z ograniczon^ amplitude taki 
sygnal, ktory powoduje maksymalne znieksztalcenie zl . Wyrazenie (36) 
uzyskuje maksymaln^ wartosc wtedy, gdy funkcja podcalkowa b^dzie 
stale dodatnia. Wynika stad, ze znieksztalcenie (36) maksymaln^ wartosc 
uzyskuje w przypadku sygnalu, ktory ma stal^ amplitude ±Emi zmienia 
znak w zerach funkcji k(r) — ko^('r— to). Dla ukladu, w ktorym funkcja ta 
nie posiada zer dla t >0 , najbardziej niekorzystnym sygnalem, ze wszyst- 
kich sygnalow z ograniczbn^. amplitude, b^dzie sygnal o ksztalcie skoku 
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Em l(t). Rezultat powyzszy uzasadnia celowosc stosowania sygnalu jednost- 
kowego jako sygnalu wzorcowego przy badaniu aparatury elektrycznej. 


3. SYGNALY MODULOWANE 

Sygnaly elektryczne stosowane w praktyce wystQpuj^ najcz^sciej 
w postaci modulowanej ^ 

e (t) — eo (t) cos ojo t = 92e { eo (f) e , 

gdzie e„(t) okresla charakter zmian obwiedni sygnalu modulowanego, 
natomiast coscoot=%{e^"«^} — reprezentuje tak zwan^ fal^ nosn^. 

Poniewaz zwykle interesuje nas obwiednia sygnalu modulowanego 
(fala nosna zostaje wyeliminowana w detektorze) wyznaczymy obwiedni^ 
sygnalu na wyjsciu ukladu liniowego. 

Wyznaczmy najpierw reakcj^ ukladu na sygnal 

■ eo(t) e‘ '’“‘ . 

Widmo tego sygnalu wyraza si^ wzorem 

E(ia)) = E„,[i (o) — Wo)] , ^ 

gdzie 

Eo(tw) = 7{eo(t)) , 

sygnal zas na wyjsciu ukladu wyraza si^ jako 

^ +00 

— r Eq [i {co - a>o)] K (i co)e*"* dco . 

— oo 

Zatem reakcja ukladu na sygnal 

ze wzgl^du na liniowosc ukladu wyniesie 

+ 00 


u (t) = 92e j E, [i (w - Wo)] K (i w) e‘“‘da) 


(37) 


Jesli sygnal modulowany dziala na ukiad o wlasnosciach selektywnych 
(na przyklad wzmacniacz rezonansowy nastrojony na czQstotliwosc nosn^) 
wygodniej jest poslugiwac si^ poj^ciem funkcji przenoszenia, w ktorej 
cz^stotliwosc ft) odmierza siQ od cz^stotliwosci rezonansowej coq,, to jest 

K(ift))=Kob[i(ft.> — ft)o)]=Kob(i^) . (38) 

Podstawiajqc (38) do wzoru (37) otrzymujemy 

-j-oo 

— - I Bo (i^)Kob m =«o (t) cos a>^t , 


Sanitized Copy Approved for Release 2010/04/07 : CIA-RDP80T00246A033700080001-9 


Sanitized Copy Approved for Release 2010/04/07 : CIA-RDP80T00246A033700080001-9 


R. Kulikowski, A. Ry bar ski 


Arch. Elektrot. 


gdzie 


Uo (t) = — r Eo (i^) Kob (i^3) dQ . 

2t.J 


■\Vz6r (39) mozemy fowniez wyrazic za pomoc^ splotu obwiedni sygnaiu 
wejsciowego e(t) z funkcj^ przejsciow^ obwiedni 

■' + OC 

2^J iQ , 

. • , .—GO ■ ^ , - 

a mianowicie . • 

CO ■ 

Uo (t) = — I 6o (t—r) hob (t) dt . 

. dt J 

0 

Wzor ten ma analogiczn^ postac do stosowanych w rozdzialach 2 i 3 wzo- 
row na splot sygnaiu e(t) z funkcj^ przejsciow^ k(t). Wynika st^d, ze 
Noptymalne charakterystyki cz^stotliwosciowe wyznaczone metod^ podan^ 
w rozdzialach 2.3, 2.4, 2.5 i odnosz^ce siQ do ukladow bez modulacji (tzw. 
videocz^stotliwdsci) mog^ bye w prostf sp'osdb przetransformowane (przez 
translacjQ czQStotliwosci i2==co — cog) na optymalne charakterystyki ukla- 
dow z modulacj^. Przy transformaeji tej ksztalt charakterystyki hie ulega 
zmianie. 

4. FILTRY „IDEALNE“ 

I • • 

W urz^dzeniach telekomunikacyjnych z tak zwanym czQstotliwoscio- 
wym rozdzieleniem sygnalow stosuje siq filtry, ktore przepuszczaj^ tylko 
okreslone (przez normy) pasmo czi^stotliwosci i obcinajq czQstotliwosci 
iez^ce poza pasmem przepuszczanym Przy projektowaniu takich filtrow 
stykamy si^ z koniecznosci^ wyznaczenia charakterystyk zapewniaj^cych 
minimalne znieksztalcenia sygnalow przy ograniezonym' wplywie zaklo- 
ceh. 

Obecnie wyznaczymy optymalne charakterystyki fiitru idealnego za- 
pewniaj^cte przeniesienie sygnaiu wzorcowego l(t) z minimalnymi znie- 
ksztalceniami 

+00 . 

>= r[hob(t)-i(t-to)Nt (40)' 


- Scisle bibr^c, uklady te hie mog^ bye zrealizoM^ne, gdyz" charakterystyki ich 
nie spelniaj^ na ogol warunku k(f)==0 dla t<0, Jednakze charakterystyki ukladow 
selektywhych mbg^ bye z wystarczaj^c^ dla praktyki dokladnosci^ aproksymow’ane 
za pomoc^ charakterystyk daj^cych sic realize wac 
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przy warunku dodatkowym (9) E = const oraz przy zalozeniu, ze widmo 
sygnalu wyjsciowegp. znika przy \Q\> Q^. 

We wzorze (40) za dolny przedzial caikowania przyjplisrhy — oo, gdyz 
sygnal hob (t) dla filtru idealnepo moze nie znikac przy t<0. Poniewaz 
widmo sygnaiu wyjsciowego 


Kob(i/3) 1 


/lob dt, 


w mysl zalozenia rowna siQ zeru na zewn^trz pasma cz^stotliwdsci prze- 
puszczanych przez filtr idealny, mozemy. napisac 

+ 0O ■ , , ‘ ■ 

— rhob(t)e-‘‘'‘dt=0, dla \Q\>Q^, (41) 

2'kJ 

—OO 

przy czym Qq jest gorn^ cz^stqtliwosci^ pasma przepuszczanego: 

Zagadnienie polega zatem na wyznaczeniu funkcji h{t ) , ktora zapewnia 
minimum funkcjonatu 

przy czym na mozliwe rozwi^zanie naklada siq dodatkowy warunek (41). 

Zagadnienie powyzsze jest zagadnieniem wariacyjnym na ekstremum 
warunkpwe. Zagadnienie takie mozna rozwiqzywac metod^ Lagrange’a, 
zwan^ tez metod^ nieokreslonego mnoznika funkcyjnego. 

Zgodnie z powyzsz^ metod^, pomnozmy wyrazenie (41) przez nieokre- 
slon^ na razie funkcjQ i scalkujmy otrzyman^ funkcji w przedziale, 

ktory le'zy na^^zewn^trz pasma czqstotliwosci przepuszczanych. 

Otrzymamy wtedy zamiast (41) warunek nastqpuj^cy 


+OC 


oo +00 




(t)e-‘«‘d£ d.O=o/ |12|;>4- 


Zmieniajqc w.e wzorze powyzszym kolejnosc caikowania otrzymujemy 

'1-00 oo 

— 00 — 00 Of, . 

W mysl metody Lagrange’a znajdujemy teraz ekstremum bezwzgl^dne 
funkcjonalu 

H = Hi+d‘=I+D+XiE. (42) 

Przyrowriuj^c do zera pierwsze wariacje wzgl^dem hob(t) oraz A (13) wy- 
razenia (42), otrzymujemy nastppuj^ce rownania: 
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iij 


-wt ^ 


“ ^ob i ^ dO + \ X(Q)e * dQ' 


-J A(i3)e'‘''‘d/3j=0, 


(43a) 


2^. 


hob(t) = 0 , dla \Q\>Qq, 


(43b) 


Poddaj^c rownanie (43a) transformacji Fouriera oraz uwzgl^dniaj^c •(43b) 
otrzymujemy 


luy \ iQ 


m. 


dla I 1 < /3o • 


Stgd oraz z (43b) otrzymujemy wreszcie-; 


Kob{iO) = 


-iOt 


— Aj (iQ)^ + (0) + 1 

0 -, 


\b\<Q,, 


|i3|>/3o. 


We wzorze tym (0) moze bye okreslone na przyklad z rownosci 


4-^0 

ii;,(o)= jKobmdQ, 


Aj zas z warunku izoperymetryeznego (9). 

Dla Aj = 0 otrzymujemy charakterystyk^ [2] 


Kob(iP^= 



.e-‘"‘ dla , :|i3|<i2o, « 

0 dla \Q\>Qo. 

Charakterystyce tej odpowiada prostokatna charakterystyka amplitudy 

i_ |1 dla \Q\<Qq, 

lo dla \Q\>Qa-, ■ 


A(i3) = |Kob(ifi)|=jJ 


oraz liniowa charakterystyka fazowa. 
Charakterystyce tej odpowiada 


hob (t) 1 Si Qq (t — to) , 

2 TZ 


gdzie Si oznaeza sinus calkowy. Charakterystyka powyzszego typu,' 
rnimo ze nalezy do charakterystyk nie daj^cych si^. zrealizowac, moze bye 
z dostateezn^ dla praktyki dokladnoscia aproksymowana przez charakte- 
rystyki ukladow mozliwych do realizaeji. 
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Zakladaj^c na przykiad, ze funkcja tlumienia ukladu wyraza siq 
wzorem 


1 


Kob(ix) tJ, fin 


(44) 

fc = l *^^0 

mozemy wyznaczyc takie rozstawienie zer x^-tej funkcji, ktore zapewnia 
najlepsze przyblizenie (w sensie Czebyszewa) charakterystyki idealnej 

=:! w przedziale (— 1,+1) oraz dostatecznie duze tiumienie 

I Koi>(ix) I 

poza pasmem przepuszczania (x>l). W tym celu zaiozmy, ze kwadrat 
modulu funkcji (44) jest proporcjonalny do 

Aq Li T2n {x) , 

gdzie T 2 n{x) — oznacza wielomian Czebyszewa 2n-tego rz^du, 

L — maksymalne odchylenie tego wielomianu od Ao= const 
Jesli nierownomiemosc tlumienia w pasmie przepuszczania ukladu 
wynosi adB (rys. la), to mamy 




e2“+ 1 


L=. 


e2a_- 1 


(rys. lb) 


Wprowadzaj^c 

tlumienia 


wspolczynnik 


’=1/— 

V e“+r 


ktory rowna si^ zeru dla a— 0 
i d^zy do je^osci p^rzy bardzo 
duzym a ; wielomian (44) mozna 
wyrazic jako ^ 


(l__T2n)2 


[l + r^”i-2T2”T2n(a:)]- 


(45) 



Hys. 1. Aproksymacja charakterystyki 
f iltru idealnego , 


Aby wielomian ten by! propor- 
cjonalny do wielomianu (44), zera 
ich powinny si? pokrywac. Zera te b^d^ lezaly w ogdlnym przypadku na 
plaszczyznie zmiennej, zespolonej x=z-{-iy . Wprowadzajqc wspolrz^dne 
eliptyczne wedlug wzorow 


» Funkcje. przenoszenia takiego typu. posiadaja na przyklad uklady zlozone ze 
wzmacniaczy rezonansowych, zawierajacych jeden lub dwa.obwody w kazdym stopniu. 
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COS (p, 


x= — (ge*’’ + — e"*’’ 
2 \ Q 


do wzoru (45) otrzymujemy * 

\i+T*^ + T^^Iq 

hT2")2 L A 


2n fyi2n<p . 


1 


-i2n(p 




(1+T2«)2 

Nie trudno teraz zauwazyc, ze funkcja powyzsza rowna siQ zeru dla 
p= — , 2nq> = {2k—l)T^ , Jc= 1 , 2 , . . . 2 ti. 

T 

wynika stqd, ze zera wielomianu (44) powinny lezec na elipsie w punktach 
o wspolrzQdnych 

1 /^ , 1 \ 2?c-l 

Zfc= — TH cos X, 

2 \ T 2n 




2k- 1 

sin^ TT. 


2n 


Przyjmuj^c zatem okreslon^ wartosc tlumienia a, mozemy ze wzorow 
powyzszych wyznaczyc optymalne polozenie biegunow Xjc funkcji prze- 
noszenia Kob(^^), a nast^pnie okreslic optymalne parametry ukladu pro- 
jektowanego. 
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* Wzor ten wynika z nast^pujacych zalezno^ci: 

= — [ (x + l/a:* - 1)^« -i(x—l/x^~ 1)^"] , 

" / 2 

1 / . 1 . \ - 

j/x^-l = Y , 

x-i-y^x^— l=Qe^^ y X— — 1= . 


\ 
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' ^ 

' P.,HynHKOBCKH, fi. PbIBflPCKH 

^ KPMTEPMn KBAAPATMHECKIIX' MCKAXCEHMH nPM OrPAHMHEHHOlt 

MOmHOCTM nOMEX ' 

P e 3 lo M e 

B KanecTBe KpMTepMH MCKaaceKWH cnmajia e(/) b JIHHei^HO^^ CMCTeiwe c nepexo^- 
ROM (J^yHKaeii h{t) npMHMMaeTCH Bbipax^eHMe I (6), r^e x(t)— jiMHewHbm chyHKUMO- 
Haji c(i), Ecjim cneKTpajibHyK) njiOTHOCTb BXO^Hbix my mob mojkho cHMTaTb B.nojioce 
nponycKaHMfl 4)MjibTpa BejinHWHOii nocTOHHHOM to cpe^HHH MomnocTb Bbixo;i;Hbix 
myMOB BbipajKaeTCH (JjopMyjioii (9). noKasaHO TaKxce, hto cpe^^w Bcex noMex c orpa- 
HMHeHHOM noBepxHOCTbK) ^ (7), sHeprHH noMex Ha Bbixo;aie rjDWJibTpa ne npeBbimaer 
BeJIMHMHbl (9). 

Bonpoc CBO^^MTCH 3 aTeM K onpeflejiCHMio MMHMMajibHoro 3HaHeHMfl MCKaxcenMii (6> 
npM ;^o6aBOHHOM ycjioBMM E= const (9), (10). TTpMpaBHiiBaH nepByio BapMaiimo <i>yHK- 
UMOHaJia W=I + ^E no h(t) k nyjiK) nojiynaeM ypaBneHMe (15), oAHoanaMHOCTb 
pemeHMM KOToporo MCCJieflOBana b rjiase 3.2. B rjiaise 3.3 npwBO^^MTCH pemeHne ypaa- 
HenwH (15), onpe^ejiHioiuee onTHMajibnyio nepexo^Hyio 4)yHKUM}0 CMCTeMbi. SaxeM 
paccMOTpeno KOHKpeTHbie npwMepbi npw bxo^hbix curnajiax 

c(t)=e-S x{i) = e(t-\) H e(i) = l(t), x(t) = l(t“-to) . 

B rjiase 4 noKasano, hto nojiyneHHbie peayjibxaTbi mojkho nepenecTn npw; 
noMOiUM TpancJiHUMM HacTOTbi Ha cjiynan aMnjiMTyflHO-rMOAyJiMpoBaHHbix cwrHajios. 

SaxeM onpe^ejiCHO xapaKTepncTiiKM cMcxeMbi, 'odecneHWBaiomeM MWHMMajibHbie 
MCKajKeHMH eAHMHHoro CHrnajia npw E = const m floSaBOHHOM ycjioBMM, KOTopoe- 

CBOflMTCH K TpeOOBaHMK), HTOGbl HaCTOTHbie xapaKTCpMCTHKM CMCTeMbi MCHCSaJIM BHC 
nojiocbi nponycKaHMH. B cjiynae jE=0 onTMMaJibHaa xapaKxepMCTMKa aMnjiMTyinbi: 
MMeex cJ)opMy npaMoyrojibHMKa, kotopbim. ksk sto noKaaano, moxcho annpOKCMMMpo- 
Baxb npM noMOiuM (t)yHKqMw saxyxaHMH (44) c jiioSom TOHHOCTbio. 

R. KULIKOWSKI A. RYBARSKI 

CRITERION FOR ROOT SQUARE DISTORTION AT LIMITED NOISE POWER 

Summary 

As a criterion for distortion of the signal e(t) in a linear sj^stem having a transfer- 
function h(t), we take expression 1(6), where x{t) is a linear operation of the func- 
tion c(t). If the spectral density of the input noise power to have a, constant 
value ^ in the frequency interval included in the filter characteristic, then the 
average output noise power will be equal to (9). It is also shown that among all dis- 
turbances with a limited area y. . (7), the output noise energy does not exceed the 
value given by formula (9). 
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The problem involves determining the minimum value of the distortion (6) at 
the subsidiary condition that E=const (9), (10). To first variation of W=I+AE with 
respect to u(t) being set equal to zero, we obtain Eq. (15). The question of wheter 
the solution of Eq. (15) is single-valued, is considered under point 3.2. The solution 
of Eq. (15) defining the optimum transfer function of the system is given under 
point 3.3. Next, specific examples with input signals c(t) = e~S x{t) = e(t-t^) and 
€(t) = l(t), a:(t) =l(t— fo) are discussed. 

» In section 4, it is shown that the results presented in the previous sections can 
be applied to amplitude modulated signals, by a translation of frequency. 

At last, the authors determine the characteristics of a system which distorts, to 
a minimum extent, a unit signal with E= const' and whose frequency characteristic 
vanishes outside a defined frequency band. In the case o,f .E=0, the optimum ampli- 
tude characteristic has a rectangular shape. It is ‘ shown that tthis shape can be ap - 
proximated to any desired accuracy with the help of the damping function (44). 
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